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Abstract 

A fourth order generalized f(R) gravity theory (FOG) is considered 
with the Einstein-Hilbert action R+aR 2 +bR IMU R^ l/ , R^ v being Ricci 's 
tensor and R the curvature scalar. The field equations are applied to 
spherical bodies where Newtonian gravity is a good approximation. 
The result is that for < a ~ — b « R 2 , R being the body radius, 
the gravitational field outside the body contains two Yukawas, one 
attractive and the other one repulsive, in addition to the Newtonian 
term. For a ~ —b » R? the gravitational field near the body is zero 
but at distances greater than \/a ~ \f— b the field is practically New- 
tonian. From the comparison with laboratory experiments I conclude 
that yfa and y/— b should be smaller than a few millimeters, which 
excludes any relevant effect of FOG on stars, galaxies or cosmology. 

1 Fourth order generalized f(R) gravity and 
quantum vacuum 

In recent years a great effort has been devoted to extended gravity theories 
which modify general relativity. The main motivation was the search for 
physical explanations to the observed accelerated expansion of the universe 
and other astrophysical observations, like the flat rotation curves in galaxies 
PP . It is plausible to derive the extended theory from a generalized Einstein- 
Hilbert action 

S =l- f d A xyf^j (R + F) + S maU (1) 



2k 



1 



where F should be a function of the scalars which may be obtained by com- 
bining the Riemann tensor, R^ v \ a) and its derivatives, with the metric tensor, 
g^ v . The theory derived from the particular choice F(R), where R is the Ricci 
scalar, has been extensively explored under the name of f (R)- gravity [2\. But 
it is possible to consider more general forms of F, e. g. depending on the 
Ricci tensor in addition to the curvature scalar. The most simple exten- 
sion of general relativity seems to be fourth order gravity (FOG), where the 
functional F has the form 

F = aR 2 + bR^R^, (2) 

where G^ u is the Einstein tensor, R^ v the Ricci tensor, R the curvature scalar 
and k is 8ir times Newton's constant, a and b being two constant parameters 
with dimensions of length squared[3], [3J. 

Another approach to FOG comes from the assumption that the energy 
of the quantum vacuum is not zero in curved spacetime[5], [6]. This fact 
may be taken into account by adding a new stress-energy tensor, T™ c , to the 
matter one, T™ at , so that Einstein's equation becomes 

G,u = R, v - \g»»R = k (r™ at + r;r) , (3) 

Eq.((3]) may be derived from the action (JTJ , the tensor T™ c coming from the 
functional F. In not too strong gravitational fields a plausible form of F is 
given by eq. fl2]) [3] . In summary, we see that FOG may be seen as either an 
extension of general relativity or a quantum vacuum effect. In the latter case 
the tensor T™ c appears on the right side of eq.((3]) , in the former it would 
appear on the left. But in both cases it gives rise to the same physical theory. 
In this paper I will use a language corresponding to the quantum vacuum 
assumption. 

The modifications derived from F should be small in weak gravitational 
fields, where GR is valid, but might be relevant in more strong fields like 
those existing in compact stars or the early universe. In order to study that 
possibility we should find the range or the parameters a and b compatible 
with known data, which is the purpose of this paper. 

2 Field equations 

The tensor field equation derived from the functional fl2]) may be taken from 
the literature [4J. I shall write it in terms of the Einstein tensor, G^ u , rather 
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than the Ricci tensor, R^, and in a form that looks like the standard Einstein 
equation of general relativity eq.flSJ). That is 

/~i _ h,n-\ rp _ rpmat i rpvac 



kT™ c = -{2a + b) [V^V.G - g^UG] - 2 (a + b) 



-b 



2G°G. 



nr - -g^uG\ a G a — Vo-V^G^ — Vo-V^G^ + OG^ v 



(4) 



with obvious notation. 

We are interested in static problems of spherical symmetry and will use 
the standard metric 



exp ((3 (r)) dt 2 + exp (a (r)) dr 2 + r 2 dVt 2 . 



(5) 



Thus G^ v (r) and T™ ai (r) have 3 independent components each, so that 



including a (r) and f3 (r) there are 8 unknown functions of r. On the other 
hand there are 8 equations, namely 3 eqs.(|l|) , 3 more equations giving the 
independent components of G^ v in terms of a and (3 and 2 equations of 
state relating the 3 independent components of T™f*. I shall assume local 
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isotropy for matter, so that one of the latter will be the equality 
(= T^ at in spherical symmetry.) In principle the remaining 7 coupled non- 
linear equations may be solved exactly by numerical methods. 

Before proceeding, a note about the signs convention is in order. As is well 
known different authors use different signs in the definition of the relevant 
quantities. Here I shall make a choice which essentially agrees with the one 
of Ref . [2] . It may be summarized as follows 



5-00 = -exp/3,G 



R^v ^gpLvR 



hT T° 
hjl iivi 1 o 



-p. 



After that I shall write the three independent components of eq.( 
notation 



mat)o 



-p, T 1 = p, Tl = q,T» = T=p + 2q-p, 

'Pmati (^mat)i = (^mai)2 = K^rnatJs = Pmat- 



(6) 

using the 
(7) 



In the following I will name p, p and q the total density, radial pressure and 
transverse pressure respectively, whilst p mat and p mat will be named matter 
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density and pressure respectively (remember that we assume local isotropy 
for matter, that is equality of radial and transverse matter pressures.) The 
differences p — p mat ,P — Pmat and q — p mat will be named vacuum density, 
radial pressure and transverse pressure, respectively. 

After some algebra I get for the components of the tensor eq. (j4j) as follows 



" d 2 T 

Pmat = -p+(2a + b)e- a —Tj- 

I dr z 

+ (a + b)k(^T 2 + 2Tp) + b 



2 1 , 

a 

r 2 



dT 

dr 



1 



-Ap + 2kp 2 - -k [p 2 +p 2 + 2q 2 ] 



+oexp(— a) 



-^-(p + q)+(lot'f3'-P")(p + p) 



(8) 



Pmat = p-(2a + b)e- a [^ + ^A^ + (a + b)k(±T 2 -2Tp) 



Ap + 2kp 2 - -k [p 2 +p 2 + 2q 2 ] 



+b 

+bexp(— a) 



(9) 



d 2 T 
dr 2 



Pmat = q- [2a + b)e 

+ (a + b)k(^T 2 - 2Tq) + b 



+&exp(— a) 



1 1 Ql 1 A dT 

- + -B a' — 

r 2 2 dr 



Aq + 2kq 2 - -k [p 2 + p 2 + 2q 2 ] 



V a' 




\ r 


r 2 ) 



r 



Addition of these 3 equations gives the trace equation, that is 

Tmat = 3p mat - p mat = T - (6a + 26) AT, 
where A is the Laplacean operator in curved space-time, that is 



A = exp(— a) 



d 2 



dr 2 



1 / 



+ - + ^' - ^ - 



d 



dr 



(10) 



(11) 



(12) 
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The quantities are related to the metric coefficients a and f3 and their 
derivatives, hence to p,p and q, that is 

2m a' m — 4irpr 3 . m + 4irr 3 p 

exp(-a) = 1 >1T = — — o =2 ~ — o ' 

r 2 r z — 2mr r z — 2mr 



„ 87rr 2 (rp + rp + 3p') 4 (m + 47rr 3 j9) (r — m — 47rr 3 p) 
r 2 - 2mr (r 2 - 2mr) 2 



where I have used units = 8tt, c = 1 and the radial derivative of a (/?') is 
labelled a' (/?"). The mass parameter m is defined by 

m = Anx 2 p(x)dx. (14) 

The condition that Einstein tensor, G^ u , is divergence free leads to the hy- 
drostatic equilibrium equation, that is 

|.fc£)_^ (p+p) . (15) 



3 Terrestrial constraints on fourth order grav- 
ity (FOG) 

The theory derived from F = in the action (CQ) , that is general relativ- 
ity, is known to give good agreement with observations for a wide range of 
intensities of the gravitational field (that is curvature of spacetime.) As a 
consequence the corrections due to finite, nonzero, values of the parameters 
a and b should be below the uncertainties in the data in that domain. In 
a previous paper [3] I considered the problem, but there are additional con- 
straints not taken into account there, which makes necessary a more detailed 
study 

In order to study the constraints on a and b derived from terrestrial and 
solar system observations I shall start solving eqs.(jHJ) to (|15p for spherical 
bodies - like the Sun, the Earth or a laboratory sphere of metal - where 
Newtonian gravity is a fairly good approximation. Corrections to Newtonian 
gravity coming from FOG have been obtained via the Newtonian approxima- 
tion of the field eqs.® to (|T0l) [4j. [?]. Here I shall calculate the gravitational 
field by solving directly the field equations, with appropriate approximations. 
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For the Earth the parameters a and /3 of the metric (??) are very close 
to unity and terms like GG V or g^G 2 are smaller than the main term, G v , by 
about 

ak Pmat/c 2 ~ bk Pmat / C 2 < 10" 26 , 

for a ~ b < 1 m 2 . In addition the matter pressure, p ma <, is negligible in 
comparison with matter density, p mat . As a consequence eqs. (JSJ) to ffTUj) may 
be approximated by the following 



d 2 p 2 dp 
<ir 2 r dr 



+ (2a + 6) 



d 2 T IdT 

dr 2 r dr 



T~ P Pmati 



-b 



d 2 p 2 dp 4 



dr 2 
2dg 
r dr 



+ 



r dr 
2 



+ (9 - P) 



, 2 c?T 
+ 2a + 6 — —-p 
r dr 



(16) 
(17) 



(2a + 6) 



d 2 T IdT 

dr 2 r dr 



g = 0. (18) 



Eqs. f|T6|) to ffTB"]) toghether with eqs.f fT3|) to f|T5l) , plus the equation of state, 
are a system of coupled differential equations whose solution is involved. 
However a great simplification is possible if we assume that the corrections 
due to finite values of the parameters a and b would modify but slightly 
the function p mat (r) with respect to the results obtained from a Newtonian 
treatment. Thus we may take p mat (r) as given, which decouples eqs . (ITS]) to 
(fl8"|) from the remaining ones. Still the three eqs. (ITE]) to (jTSl are coupled 
amongst themselves, but from them it is possible to get two decoupled ones. 
In fact if I add eq.(!T6l) , eq. (fT7l) plus twice eq.(TT8"l) I get the trace eq.lJTTT). 
which may be rewritten 



r d 2 T 2 dT 

dr 2 r dr 



(6a + 26) 



— T 



-T, 



mat 



P 



mat' 



Subtracting this minus three times eq.( TT6|) I obtain 
"d 2 (T + 3p) 2d(T + 3p) 



+ T + 3p = 2p 



mat ' 



dr 2 r dr 

Now the general solution of the trace eq.( JT9l) is trivial and I get 

1 



(19) 



(20) 



T(r) 



y/Qa + 2br 



A 

H — exp 

r 



sinh 



VQa + 2b 



r — z 

y/6a + 2b 
B 

H exp 



ZPmat 



\z] dz 



V6a + 2b 



(21) 
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The integration constants A and B may be got from the boundary conditions, 
that is the function p (r) should be finite at the origin and go to zero at 
infinity, which gives 



1 



R 



A = —B = / exp zPmat ( z ) dz, (22) 

where R is the radius of the body, that is, the radius of the matter distri- 
bution. Note that a part of the total mass of the body lies in the region 
without matter, associated to the vacuum density. Eq.( )20|) may be solved 
by a method similar to the one used for eq.f ll9p , which gives the function 
p (r) . Once we know the function T(r) the solutions of eqs.(TT7l) and (fT5|) are 
straightforward. 

I am interested in the functions T (r) , p (r) , p(r) and q{r) outside the 
body (i.e. r > R, where p mat = 0) and for the particular case where p mat is a 
constant inside the body ( which is a good approximation for the Earth or a 
metallic sphere). Thus the z integrals in eqs.f l2"T|) and (|2"2"1) may be performed 
analytically, and similarly in the solution of eq. ( 12"U|) . In order to simplify the 
notation I shall introduce the following dimensionless variables 

r R r R 

x= - ,X= - ;Z/ = ^_y = ^_ (23) 

V6a + 2b VQa + 2b y^b ^b 



and the mass parameters 

— [X - 1 + (X 4 -• ' i ; - 2V , 



M x = ^[X-l + (X + l) exp (-2X)} , My = ^ [Y - 1 4 (Y 4 1) exp (-2Y)} 



(24) 



Thus I get for r > R 



r M = - I^T26)7 exp < x -*>' ' 25 > 

C( r ) = TT7 — K 6a + 2b Y 1 M * exp (X - i) + 2 |6P 1 M y exp (y - I/)] . (26) 
IZnr 

Hence it follows 

p + 2q = T + p=-^— \\b\~ 1 M v exp (Y — y) — (6a 4 2b)~ l M x exp (X - a;)] . 
6vTr 

(27) 
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In order to get p and q separately we shall solve the equation resulting from 
the subtraction of eq. (fTSl) minus eq. CTTj) in the region r > R, that is 



|6| 



d 2 h + 2dh 6h 
dr 2 r dr r 2 



lW lt£T d 2 T, . . 

~h = 2a + 6 -— - — - (28) 
r dr dr z 

(2a + 6) M x ( 3 3 \ 

' 1 + - + -» exp {X - x) , 



47r(6a + 26) 5 /2 y x x 2 
where eq. (l25|) has been taken into account. Now I use the ansatz 

( 3 3 \ / 3 3 \ 

h = A I 1 + - + — 1 exp (X - x) + B I 1 -\ h ) exp (Y - y) , 

\ x x J \ y y J 

which inserted in the left side of eq. (|28|) gives 
This leads to 

A( \ b \ _ 1)= ga + b) M x 
{ Qa + 2b ' 6a + 26 12tt' 

whence I get 



h = -^(l + - + ^)exp(X-x)+B(l + - + - 2 )e W (Y-y). 

±27T \ x x z j \ y y 



Combining this with eq. (l27j) I obtain 
P = g^: [\b\- 1 M y exp(Y-y)-(6a + 2b)- 1 M x exp(X-xj\ 

+ ir( 1 + - + ^) ^v{X-x)-2B (l+* + ^)e W (Y-y). 

d7t \ x x j \ y y J 

The integration constant B is fixed by the condition that eq. (1151) is fulfilled 
with the approximation of neglecting the last term. That is 

dp 2(q-p) d , 3 v 2 

— = — (r 3 p) = r z (2q + p) . 

dr r dr 

Finally this leads to 
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P( r ) = 7T^i x + 1] ex P {X-x)-^[y + V\ exp (Y — y) . (29) 



All results obtained up to now for the solution of the differential eqs. tjlb]) to 
( 12 Op have assumed that both 6a + 2b and —b are positive. If one or both 
quantities were negative we would obtain sinus or cosinus functions rather 
than exponentials, which would clearly violate empirical facts. Thus the 
parameters are constrained to the range 

< -b < 3a. (30) 

The remaining bounds are derived in the following. 

The vacuum density eq.( l26|) gives rise to a mass distribution which I 
calculate as follows. Firstly I show that the total mass of the body is the 
same as the mass of matter, that is the total vacuum mass is zero. This 
follows trivially if we integrate the two sides of eq. fll6p after multiplication 
times the volume element. In fact we obtain 

f°° 2 2 , f°° d 2 . . , r°° d , s , 

/ V pAnr dr = An / r ~r^ vP) " r = ~^ 7r / ~r vP) " r = 0> 
Jo Jo dr Jq dr 

where I have performed an integration by parts and taken into account that 
rp —7- for r — >■ oo. A similar result holds true in the integral of V 2 T. Thus 
the mass asociated to the total density p equals the mass asociated to the 
matter density p mat - In our language we may say that the total mass of the 
quantum vacuum f associated to T™ c , see eq.flU)) is zero. Nevertheless the 
vacuum density changes the spatial distribution of the total mass. In fact, we 
may obtain the (vacuum) mass distribution outside the body by integrating 
eq.tES]) between R and r. I get 

M ext (r) = An r 2 pdr 
Jr 

= — - [X + 1 - (x + 1) exp (X - x)} + ^ [Y + 1 - (y + 1) exp (Y 

Hence it follows that the total mass enclosed in a sphere of radius r is (see 
eqs.((23D) 

M (r) = M - — - (x + 1) exp (X - x) - ^— ^ (y + 1) exp (Y - y) , (31) 
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where the condition M (oo) = M has been taken into account. Remember 
that M (oo) — M corresponds to the vacuum mass, which is zero as shown 
above. The effect is dramatic for a spherical body with a radius R << 
y/a, y/\b\. In this case X, Y « 1, M x ~ M y ~ M which implies M (r) ~ 0, 
that is the total mass in the interior of the body is zero. That is the vacuum 
mass (negative) cancels the matter mass (positive) in the interior of the body. 

Now we may calculate the gravitational field, g, near the Earth surface, 
which we should identify with —1/2 times the quantity defined in (fi~3l) , 
where we may neglect m « r in the denominator. Thus from eqs.( !26]) and 
(|29|) I get, to lowest order in y/6a + 2b and yffe[, 



Q 

— - [M(r) + 47rr 3 p(r)] 



4GM V 

H ^ 

3r 2 



(y + l)exp(y- j/192) 



where G is Newton constant. We see that when \/6a + 26, « R eq.(l32l) 
corresponds to Newtonian gravity plus a small correction consisting of two 
Yukawa-type terms, one of them attractive, the other one repulsive. Although 
these results have been derived for bodies with constant density they are valid 
for any celestial body because the parameters M x and M y depend only on 
a small region near the body's surface where the density is effectively a 
constant. 

In sharp contrast for a small sphere where \/6a + 26, y/\b\ » R and 



M x ~ 



My ~ M (see eqs.((24])) we obtain 



CM 



exp (X — x) — — exp (Y — y) 



3r 2 



r-i?) 



V6A + 26 



<< 



GM 



for r not much larger than both \/6a + 26 and ^/|6[. Now performed labo- 
ratory experiments have shown that the field g agrees fairly well with New- 
tonian predictions for bodies greater than a few millimeters [7]. Thus our 
calculation shows that these experiments exclude values of the parameters a 
and |6| greater than a fraction of squared centimeter. 



10 



4 Conclusions 



I conclude that y/a and a/— b should be smaller than one centimeter, which 
probably excludes any relevant effect of FOG on stars, galaxies or cosmology. 
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